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Abstract 
> 

We show that a hadron gas model with continuous particle emission instead 

of freeze out may solve some of the problems (high values of the freeze out 

ON ■ density and specific net charge) that one encounters in the latter case when 

studying strange particle ratios such as those by WA85. This underlines the 

necessity to understand better particle emission in hydrodynamics to be able 



c 



to analyze data. It also re-opens the possibility of a quark-hadron transition 



^ ■ occur ing with phase equilibrium instead of explosively. 

53' 



1. Introduction 

An enhancement of strangeness production in relativistic nuclear collisions (compared 
to e.g. proton-proton collisions at the same energy) is a possible signature [|1|] of the much 
sought-after quark-gluon plasma. It is therefore particularly interesting that current data 
at AGS (Alternating Gradient Synchroton) and SPS (Super Proton Synchrotron) energies 
do show an increase in strangeness production (see e.g. [0]). At SPS energies, this increase 
seems to imply that something new is happening: in microscopical models, one has to 
postulate some previously unseen reaction mechanism (color rope formation in the RQMD 



code ||, multiquark clusters in the Venus code [§]], etc) while hydrodynamical models have 
their own problems (be it those with a rapidly hadronizing plasma || or those with an 
equilibrated hadronic phase, preceded or not by a plasma phase). In this paper, we examine 
the shortcomings of the latter class of hydrodynamical models and suggest that they might 
be due to a too rough description of particle emission^]. 

To be more precise, let us assume that a hadronic fireball (region filled with a hadron 
gas, or HG, in local thermal and chemical equilibrium) is formed in heavy ion collisions 
at SPS energies and that particles are emitted from it at freeze out (i.e. when they stop 
interacting due to matter dilution). One then runs into (at least) three kinds of problems 
when discussing strange particle ratios. 

First, the temperature (Tf out ~ 200 MeV) and baryonic potential (fib f. out ~ few 100 
MeV) needed at freeze out [|6|-[T0|| to reproduce WA85 [|ll|] and NA35 |L2| strangeness data 
actually correspond to high particle densities: this is inconsistent with the very notion of 
freeze out []. 

Second, to reproduce strange particle ratios, it turns out that the strange quark potential 
/i s must be small and the strangeness saturation factor 7 S of order 1 (this quantity, with 
value usually between and 1, measures how far from chemical equilibrium the strange 
particles are, 1 corresponds to full chemical equilibrium of the strange particles). Both facts 



x The main problem for the former class of hydrodynamical models is the difficulty to yield enough 
entropy after hadronization. 

2 While WA85 and NA35 data for strange particle ratios are comparable and lead to high T's and 
/Vs, NA36 data are different and lead to lower T's (Phys. Lett. B327 (1994) 433) for similar 
targets but a somewhat different kinematic window. However the rapidity distribution for A's 
(E.G. Judd et al. Nucl. Phys. A590 (1995) 291c) as well as A's and K°'s (J.Eschke et al. Heavy 
Ion Phys.4 (1996) 105) for NA36 are quite below that of NA35; NA44 midrapidity data for AT ± 
agree with that of NA35. 



are expected in a quark-gluon plasma hadronizing suddenly, not normally in a hadronic 
fireball 00. 

Third, using the values at freeze out of the temperature, baryonic potential and satura- 
tion factor extracted to reproduce WA85 strange particle ratios, one can predict the value 
of another quantity, the specific net charge (ratio of the net charge multiplicity to the total 
charge multiplicity). This quantity has been measured not by WA85, but in experimental 
conditions similar to that of WA85 by EMU05 [ 15| . It turns out that the predicted value is 



too high, (while it might be smaller if a quark-gluon plasma fireball had been formed) fll6lj5l. 
In what follows, we study how problems 1 and 3 are related to the mechanism for particle 
emission normally used, freeze out, and suggest that the use of continuous emission instead 
of freeze out, might shed some light on these questions. (We also re-discuss problem 2.) 
This underlines the necessity to understand better particle emission in hydrodynamics and 
re-opens perspectives (see conclusion) for scenarios of the quark-hadron transition. 

2. Fluid behavior and particle spectra 

First let us see in more details what are the two particle emission mechanisms just 
mentioned. In the usual freeze-out scenario, hadrons are kept in thermal equilibrium until 
some decoupling criterion has become satisfied (then they free-stream toward the detectors). 
For example, in the papers mentioned above where experimental strange particle ratios are 
reproduced, the freeze out criterion is that a certain temperature and baryonic potential have 
been reached. The formula for the emitted particle spectra used normally is the Cooper- 
Frye formula ||17|| . In the particular case of a gas expanding longitudinally only in a boost 



invariant way, freezing out at some fixed temperature and chemical potential, the Cooper- 
Frye formula reads 

-; -j — = -^—Tf.outm± ^2(T) n+1 exp (n/j,f. out /Tf, out )K 1 (nm±/Tf. out ). (1) 

(The plus sign corresponds to bosons and minus, to fermions.) It depends only on the 
conditions at freeze out: Tf out and /i/. D «t = Hbf.outB + Us f.outS, with B and S the baryon 



number and strangeness of the hadron species considered, and ps f.<wt(lJ>b f.out, Tf, out ) obtained 
by imposing strangeness neutrality. So the experimental spectra of particles teach us in that 
case only what the conditions were at freeze out. 



In the continuous emission scenario developed in ||lq , [T9|j , the basic idea is the following. 
Due to the finite dimensions and lifetime of the fluid, a particle at space-time point x, 
has some chance V to have already made its last collision. In that case, it will fly freely 
towards the detector, carrying with it memory of what the conditions were in the fluid at 
x. Therefore the spectrum of emitted particles contains an integral over all space and time, 
counting particles where they last interacted. In other words, the experimental spectra 
will give us in principle information about the whole fluid history, not just the freeze out 
conditions. For the case of a fluid expanding longitudinally only in a boost invariant way 
with continuous particle emission, the formula that parallels (1) is 

dN 2g r m± cosh rirppdp + p± cos (ppprdr , . 

dcpdrj — — — (2) 



dyp±dp± (27r) 2 Jv=o.5 exp((mj_ cosh 77 — p)/T) ± 1 

where 7>(p, 0, 77; v± (resp. Pf(t, 4>,tj; v±)) is the time (resp. radius) where the probabil- 
ity to escape without collision V = 0.5 is reached. V is given by a Glauber formula, 
exp[— J av re in(r')dT'} and depends in particular on location and direction of motion. We 
are using both (1) and (2) in the following. Clearly, in (2), various T and p = pbB + psS 
appear (again p s (pb,T) is obtained from strangeness neutrality), reflecting the whole fluid 
history, not just T f , out and p b f. ou t- 

So to predict particle spectra, in the case of continuous emission, we need to know the 
fluid history. To get it, we fix some initial conditions T(r ,p) = T and /if,(r ,p) = fibo and 
solve the equations of conservation of momentum-energy and baryon number for a mixture 
of free and interacting particles, using the equation of state of a resonance gas (including 
the 207 known lowest mass particles) and imposing strangeness neutrality. As a result we 
get T(r,p),fj,b(r,p) and we can use these as input in the formula for the particle spectra 
(2). The procedure is similar to that of a massless pion gas [P|P|] but is numerically more 
involved. 



An important result [jT8|JT9|] for the following is that for heavy particles, the spectrum 
(2) is dominated by the initial conditions, precisely a formula similar to (1) with freeze out 
quantities replaced by initial conditions could be used as an approximation, (particularly 
at high p±); for light particles the whole fluid history matters. To understand this fact, 
one can consider equation 2 and compare particles emitted at T(r, p)=200 and 100 MeV. 
For particles with mass of 1 GeV, the exponential term gives a thermal suppression above 
one hundred between these two tempratures. The multiplicative factors in front of the 
exponential are in principle larger at the lower temperature but do not compensate for such 
a big decrease. This is why heavy particles are abundantly emitted at high temperatures. 
On the other side for pions, the thermal suppression is only a factor of 2. This is why light 
particles are emitted significantly in a larger interval of temperatures. 

Note that since heavy particle and high p± particle spectra are sensitive mostly to the 
initial values of T and //&, the exact fluid expansion does not matter very much for them; in 
particular the assumption of boost invariance should play no part in the forthcoming analysis 
of strange high p± particle ratios]]. It would be however interesting to include continuous 
emission in, e.g. a hydrodynamical code, to obtain the fluid evolution and study pion data 
and low p± data. 

3. Particle ratios 

Once the spectra have been obtained, they can be integrated to get particle numbers, 
taking into account eventual experimental cutoffs and correcting for resonance decays. Since 
we had to specify the initial conditions to solve the conservation equations and use this 
solution as input into (2), the particle numbers depend on T , /j,bo- In contrast, for the freeze 



3 Note also that the data considered below are in a small rapidity window, near midrapidty. Were 
it not for this fact, boost invariance should not be assumed, because the rapidity distributions do 
not have this symmetry. 



out case, particle numbers depend on the conditions at freeze out, Tf_ out and fibf.out- 

We look for regions in the T , fij,o space which permit to reproduce the latest WA85 
experimental data on strange baryons |jTTj| for 2.3 < y < 2.8 and 1.0 < p± < 3.0 GeV: 
A/A = 0.20 ± 0.01, S^/S" = 0.41 ± 0.05 and S~/A = 0.09 ± 0.01 (S^/A = 0.20 ± 0.03 
follows). In fact there is no set of initial conditions which permits to reproduce all the above 
ratios. A similar situation occurs with freeze out, as noted in ||20|| . (Note that with the older 
experimental data pl| , a region in the Tf out — Hbfout space permiting to reproduce all the 
above ratios existed |6|j8]-|T0|.) 

In the comparison of our model with WA85 data we have assumed however complete 
chemical equilibrium for strangeness production. As already mentioned in the introduction, 
this is not expected for a HG. In order to account for incomplete strangeness equilibration, 
we introduce the additional strangeness saturation parameter 7,, by making the substitution 
exp usS — > jf^expusS in the (Boltzmann) distribution functions [^]. In our case, a 
priori, 7,, depends on the space-time location x, however since as already mentioned, the 
initial conditions dominate in the shape and normalization of the spectra of heavy particles 
(particularly at high mi), we take 

dN 7i S| (r„)^%- (3) 



dyp±dp ± s dyp±dp A 

with dN eq /dyp^dp± given by (2). In figure la, we see that for 7 s (r ) = 0.58, there exists an 
overlap region in the T ,/if,o plane where all the above ratios are reproduced. For the freeze 
out case, a similar situation occurs as noted in p0f , namely there exists an overlap region 
for 7 S = 0.7. 

In the freeze out case, the values of the parameters in the overlap region correspond to 
high particle densities and so it is hard to understand how particles have ceased to interact: 
this is the problem 1 mentioned in the introduction. In the continuous emission case, T and 
/ifeo in the overlap region lead to high initial densities, but this is of course quite reasonable 
since these are values when the HG started its hydrodynamical expansion. Note also that 
the value we used for < av re i > comes from using a Breit-Wigner formula for the A — it cross 
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section and computing < av re i > p a >10 at various temperatures. 



The aim of picture la is to allow an easy comparison with freeze out results such as pOf . 
however it is not physically complete: so far we have neglected hadronic volume corrections. 
For freeze out, this correction cancels between numerator and denominator in baryon ratios 
so can be ignored [|II]] but for continuous emission, since we are considering the whole fluid 
history to get particle numbers (and then their ratio), it must be included. There are various 
ways to do this. (Some of the) thermodynamical quantities for pointless particles can be 
multiplied by a factor (determined from pointlike quantities), for example 1/(1 + Vtiq), 
where V is a typical hadronic volume and indicates pointlike quantities. However there 
is no consensus on whether this correction involves a net baryonic density n p3|j2^JT0|j2 



or the total particle density n = J^jrf p6| , p7| or a compromise ||. Here we modify all 



particle densities as follows n, = n®/(l + VJ2j n ?) as m [P^ , p7|1 - Taking into account hadronic 



volumes, we get the overlap region shown in figure lb, which is shifted towards smaller T's 
and /ife's but not very different from that of figure la. Given that simulations of QCD on 
a lattice indicate a quark-hadron transition for temperatures^ around 200 MeV, it seems 
more reasonable to consider initial conditions T ~ 190 MeV and ^q ~ 180 MeV, i.e. the 
bottom part of the overlap region. The precise location of the overlap region (and exact 
value of 7 S ) is sensitive to changes in the equation of state -as we have just seen- as well as 
in the cross section or cutoff V = 0.5 in eq. 2. (As a cross check, we have also remade our 



calculations using a Walecka-type equation of state [28] and found similar results to figure 



4 In these simulations, the transition temperature depends strongly on the number of flavors Nf 
(being > 200 MeV for Nf = and < 200 MeV for Nf > 1) and somehow on the method of extrac- 
tion. Since the created plasma might remain gluon dominated with quarks away from equilibrium, 
up to the transition (this is already predicted to be the case for RHIC and LHC energies), it is not 
clear what value of the transition temperature should be used but one should probably not take 
results from lattice QCD as a rigid value in our context. 



lb.) Therefore problem 1 (whether the overlap region is physically reasonable) is taken care 
of. 

To be complete, we also examined the more recent ratios obtained by WA85 |29[ (at 
midrapidity): TF/n m± ^ 23GeV = 0.57 ± 0.41, (fi~ + TT=)/(Er + ^) m± >2.3GeV = 1.7 ±0.9, 
K°JA px>1 . QGeV = 1.43±0.10, K°/A p±>1 . 0G eV = 6.45±0.61 and K+/K p±>09GeV = 1.67±0.15. 
We looked for a region in the T , /j, b0 plane where ^~/^mj_>2.3Gey * s reproduced: due to the 
large experimental error bars, this does not bring new restrictions to fig. lb. We also 
calculated our value for (Q~ + Q~)/(E~ + 'E,~) m± >2.3GeV in the overlapping region and found 
~ 0.7, in marginal agreement with the above experimental values. The three ratios involving 
kaons depend on more than just initial conditions (kaons are intermediate in mass between 
pions and lambdas so part of the fluid thermal history must be refleted in their spectra), in 
particular j s (x) ~ 7s (to) ~ est may not be a good approximation and we are still working on 
this. The above experimental ratios concern SW collisions, data with SS are not so extensive 
yet but not very different |JD| so a similar overlapping region can be found. 

The apparent temperature extracted from the experimental p± spectra for A, A, S~ and 



2~ is ~ 230 MeV [I IJ. Given that we extracted from ratios of these particles, temperatures 
To > 200 MeV, we conclude that heavy particles exhibit little transverse flow, which is 
compatible with the fact that they are emitted early during the hydrodynamical expansionF]. 

4. Specific net charge 

We now turn to 

D q = (N + -N-)/(N + + N~) (4) 



5 Note that WA85 data concern highpj_ strange particles. But in [18,19|, we showed that (a simpler 
version of) our model with continuous emission and no tranverse flow, reproduces the shape of the 
NA35 S+S p± spectra for strange particles, which extend to low p± (but are not restricted to 
midrapitidy) . 
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using the continuous emission scenario. As mentioned in the introduction, for HG models 
with freeze out the predicted D q is too high, when using values of the freeze out parameters 
that fit strangeness data, e.g. Tf MUt ~ 200 MeV, fibf.out ~ 200 MeV and j s ~ 0.7. For contin- 
uous emission, due to thermal suppression, particles heavier than the pion are approximately 
emitted at T ~ 200 MeV, fi b0 ~ 200 MeV and 7 S ~ 0.49 (fig. lb), so D q so far is similar 
to that of freeze out. However there is an additional source of particles that enters the de- 
nominator of (4), namely pions are emitted at To and then on (since they are not thermally 
suppressed). So we expect to get a lower value for D q in the continuous emission case than 
in the freeze out case. (We recall that pions are the dominant contribution in iV + + N~ .) 
This would go into the direction of solving problem 3, it is still under investigation. 

5. Conclusion 

Our present description is simplified, for example we do not include the transverse ex- 
pansion of the fluid, use similar interaction cross sections for all types of particles, etc. In 
addition, we need to look systematically at strangeness data from other collaborations as 
well as other types of data such as Bose-Einstein correlations. Nevertheless, we have seen 
that the continuous emission scenario with a HG may shed light on problems 1 and 3 (dis- 
cussed in the introduction) that a freeze out model with a HG encounters. Namely, in the 
overlap region of the parameters needed to reproduce WA85 data, the density of particles 
is high and this is consistent with the emisssion mechanism, since it is the initial density 
of the thermalized fluid. We also expect D q to be smaller for continuous emission than 
freeze out. But (problem 2) the value of the strangeness saturation parameter may be high 
for a HG, particularly at the begining of its hydrodynamical expansion. However what we 
really need to get fig. lb, is that S _ /A = 7 S S~/7AA| eg . and S~/A = 7 H H~/7AA| eg . with 
7h/7a = 0.49. We expect indeed that multistrange S~ and H~ are more far off chemical 
equilibrium than singlestrange A and A so that 7= /7a < 1- The result 7 S = 0.49 arises if 
one makes the additional hypothesis that quarks are independent degrees of freedom inside 
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the hadrons so that one has factorizations of the type j\exp fi^/T = 7 s exp2/z g /Texp /j, s /T 
and 7=exp/x=/T = 7^ exp fi q /Texp2fi s /T. Therefore problem 2 may not be so serious. 

The fact that we may cure some of the problems of the HG scenario does not mean that no 
quark-gluon plasma has been created before the HG, in fact it may open new possibilities for 
scenarios of the quark-hadron transition (e.g. an equilibrated quark-gluon plasma evolving 
into an equilibrated HG with continouous emission); in particular it may not be necessary 
to assume an explosive transition || or a deflagration/detonation scenario [J3lH33|- 

But our main conclusion is that the emission mechanism may modify profoundly our in- 
terpretation of data (for example, does the slope in transverse mass spectrum tell something 
about freeze out or initial conditions?). In turn this modifies our discussion of what potential 
problems (such as 1 and 3) are emerging. Therefore we believe it is necessary to devote more 
work to get a realistic description of particle emission in hydrodynamics, fL8| , |19| being a first 
step in that direction. We remind that the idea that particles are emitted continuously and 
not on a sharp freeze out surface is supported by microscopical simulations at AGS energies 
and SPS energies 
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Note added 

After completing this paper, we learned that G.D. Yen, M.I. Gorenstein, W. Greiner 
and S.N. Yang suggested |36| another solution to problems 1 and 3 above, in term of the 
excluded volume approach of p| , for the preliminary Au+Au (AGS) and Pb+Pb (SPS) 



data. 
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FIGURE CAPTIONS 

FIG.l Overlap region in the To — fXbo plane for WA85 data, with < av re i >= 1 fm 2 a) 
without b) with hadronic volume corrections. 
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